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ABSTRACT 

We present a thorough hterature study of the most-massive star, mmax, in several 
young star clusters in order to assess whether or not star clusters are populated 
from the stellar initial mass function (IMF) by random sampling over the mass range 
0.01 ^ m ^ 150 M0 without being constrained by the cluster mass, Moci- The data 
reveal a partition of the sample into lowest mass objects (Med ^ IO^Mq), moderate 
mass clusters (lO^M© < Mod ^ lO^M©) and rich clusters above 1O^M0. Additionally, 
there is a plateau of a constant maximal star mass (mmax ~ 25 Mq) for clusters with 
masses between 10'^ Mq and 4 • IQ^Mq. Statistical tests of this data set reveal that 
the hypothesis of random sampling from the IMF between 0.01 and 150 Mq is highly 
unlikely for star clusters more massive than IO^Mq with a probability of p w 2 • 10"'' 
for the objects with Med between IO^Mq and lO'^M© and p « 3 • 10~^ for the more 
massive star clusters. Also, the spread of mmax values at a given Afoci is smaller than 
expected from random sampling. We suggest that the basic physical process able to 
explain this dependence of stellar inventory of a star cluster on its mass may be the 
interplay between stellar feedback and the binding energy of the cluster-forming molec- 
ular cloud core. Given these results, it would follow that an integrated galactic initial 
mass function (IGIMF) sampled from such clusters would automatically be steeper in 
comparison to the IMF within individual star clusters. 

Key words: stars: formation - stars: luminosity function, mass function - galaxies: 
star clusters - galaxies: evolution ~ galaxies: stellar content - Galaxy: stellar content 



1 INTRODUCTION 

Whether or not newborn stars in star clusters are randomly drawr0 from the IMF is of utmost importance for various fields 
of stellar astrophysics. For example non-random drawing which suppresses the number of OB stars in smaller clusters would 
steepen the IMF for who le galaxies, the integrated galactic stellar initial mass function (IGIMF, IXroupa fc Weidneilbood : 
Weidner fc Kroupall2005l ). A randomly-drawn IMF on the other hand, which would be equivalent to postulating the existence 



of clusters comprised of a massive star and not much more, would not (|Elmegreenll2006 Selman fc Melnick 200S ). As the bulk 



of the galactic field star populations are probably made from dissolving star clusters ( Kroupa 19951 : Lada fc Lada 19951 . 20031 : 



Adams fc MversI 2001 : Allen et al. 2007 ). understanding the stellar distribution in galaxies presupposes knowledge of the IMF 



in star clusters. A central issue on deciding whether a star cluster can be modelled in terms of random sampling from the IMF 
or not is the existence of a non-trivial relation between the mass of the most-massive star (mmax) and the star cluster mass 
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(Meci). Thus, a statistically significant correlation mmax(Moci) would imply physical processes such as self-regulation of the 
star-formation pro cess in a cluster. T he influence of the cluster mass or density on its stellar population has been studied on 
previous occasions. iLarsonl (jigsd ) and lLarsonI l|2003l ) examined the properties of molecular clouds and the stellar populations 
found within them, finding the following empirical expression between the mass of the most-massive star, rrimax, and the 
stellar mass of an embedded cluster. Mod, 

^Larson -i r)ii4-0.45 

"^max = l-2Mecl , (1) 

whi ch is shown as a d ash- dotted line in Fig. [T] 

Elmegreen (Il983l) i nvest igated a model for the formation of bound star clusters where the luminosity of the stars chosen 



from a Miller fc Scalol (|l979h IMF overcomes the binding energy of a molecular cloud. Different star formation efficiencies 
would then determine if a cloud becomes a bound star cluster or an OB association. He also found a relation bet ween mmax 
and M cci which cannot be written as an analytical equation and is shown as a long-dashed line in Fig. [1] Later, lElmegreen 
derived a different relation when assu ming a sin g le slo pe power-law IMF, ^(m), where dN = £_(m)dm is the number of 
stars in the mass interval m, m + dm, with a lSalpeterl (|l955l ) slope and solving the following two equations. 



and 

Med = 



^{m)dm 



m^(m)dm, 



but wi thout any limit for masses of the stars, mmax* ~ ( 
with a lSalpeter slope these two equations yield in. 



Here, 



Elmegreen 



/ Med y 

U X 103/ 



100, 



(2) 



(3) 



is the minimum mass. For a single power-law IMF 



(4) 



.3 X 103, 

shown as the short-dashed line in Fig. [T] 

In their numerical calculations of star-forming molecular clouds using a smoothed particle hydrodynamics code lBonnell et al 
1 20031 . [2004 ) found a relation. 



0.39 X M-^i^, 



(5) 



which is shown as a dotted line in Fig. [T] 

In a thorough st udy of star clu s ters a nd OB associations in order to determine whether or not a fundamental upper mass 
limit for stars exists. lOev fc Clarke (2005) also calculated the expected dependence of rrimax on i\fcd if the stars are randomly 
drawn from an IMF, 
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plotted as a short- dashed-long- dashed line in Fig. [T] 

Including a fundamental upp er mass limit for s tars, rrtmax* = 15OMe0 in eqs.[2]and[3l and using the canonical multi-part 

power- law IMF (Appendix 1 Weidner fc Kroupal (|2004l ) found the relation visible as a thick-solid line i n Fig. [T] 

As evident from Fig. [1] these studies arrive at a rather large range of possible mmax-Med-relations. IWeidner fc Kroupa 



I 2006h re-investigated this question by compiling a larger number of observational results from the literature and extensive 
Monte-Carlo experiments of different sampling algorithms and found evidence that there exists a n on-trivial relation betwee n 
the mass of a star cluster an d the most-massive star in the cluster, a result in principle confirmed by Selman fc Melnick ( 20081 ) . 
But they conclude that the I Weidner fc Kroupa (|2006l ) sample is biased by a size-of-sample effect. Furthermore, in the recent 
literature several claims h ave been made against such a relation arguing instead for a pure random sampling from the IMF 



in individual star clusters ( Qev et al. 2004 



de Wit et al 



2OO5I: lElmegreenlbood: barker fc GoodwinllioOTi: ISelman fc Melnick 



20081 : iMaschberger fc ClarkejboOSl l. lde Wit et all ( 20041 . 12OO5I ') find up-to 4% of non-runaway (less than 30 km/s space motion) 



O stars in isolation with no apparent cluster around them or within their lifetime if they would have been ejected from a 
cluster with a velocity of 6 km/s - indicating they formed outside a cluster. This result would of course be irreconcilable 
wit h a relation betwee n the mass of the most-ma ssive star and the mass of its parent star clu ster as has been pointed out 
by 
in 



Parker fc Goodwi: 



Weidner fc Kroup; 




and ISelman fc Melni ck (2008). Whi le Selman fc Melnickl (l2008l) argue that the sample used 
is biased against random sampling, Parker fe'coodwinl ( 200^ find that the observed 4% of 



allegedly isolated O stars would agree with random sampling from a cluster number distribution function which scales with 



150 Mp) is beli e ved to be th e fundamental upper mass limit for stars with non-zero metallicity jWeidner fc KroupalliooS : lOev fc Clark j 
I2OO5I: lFieeJl2005l : lKoenll2006l) . 
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Figure 1. Different relations between the most-massive star and the cluster mass fro m observation s, numerical calculations and theoretical 
modelling from the literature. The dash- dotte d line marks the em pirical relation bv iLarsonl l l2003t) . th e long-dashed li ne is the lElmeereenI 
lll983h relation, the short-dashed li ne is the Elmcgree™ relation, the dotted line shows the iBonnell et al. I |2003t) relation, the 

short- dashed-lonq- dashed line is the lOev &: Clarke (2005) relation and the thick-solid line is the analytical model of Weidner &: Kroupal 
1I2OO4I) . See text for more details. 



A''^, where is the number of stars in a cluster. But it should be noted here that the de Wit et al. 1 20051 ) result is an upper 
limit for O stars formed in isolation and that more in-depth observations might reduce this sample. For example, HD 165319, 
an 09. 5 I star from the de Wit et al. I 1I2OO5I ) sample of 11 stars which are indicated there as one of "the best examples for 
isolated Ga lactic high-mass star formation " has a bow-shock front an d is therefore a star ejecte d from a star cluster, possibly 
NGC 6611 (ICvaramadze fc Bomanjboosh . Additionally, according to lSchilbach fc Roseij (I2OO8I ) further 6 of the remaining 10 
stars are at distances to star clusters only slightly larger than what they may have travelled during their expected life times. 
But the current large errors of the space motion of these stars does not allow to constrain the birth places of them. 

A non-trivial mmax-Mod-relation, and therefore whether or not the stars in star clusters are randomly sampled from the 
IMF, would also give more insight and understanding of the process of star-formation. The formation of massive stars (> 
10 M0) is still not well understood with a t least two competing theor i es (competitive accretion vs. single sta r accretion) 
having been developed (jBonnell et al.lll998l , |2004| : iBonnell fc Batelbood : iTan et ahlbood : iKrumholz et ali 12009*). An m^ax- 
Mcci-relation could imply that the bulk of the low-mass stars form first and the high-mass stars later. The combined feedback 
of the massive stars would then halt further star-formation. 



Here we will show that the observed distribution of the mass, mmax, of the most-massive star in a star cluster cannot be 
drawn randomly from the IMF for clusters more massive than 100 Mq but that there must exist a physical relation between 
w^max and the birth stellar mass of the cluster, Med (the stellar content before gas expulsion but after cessation of star 
formation) . 



2 THE DATA 

2.1 Sample construction 

In order to construct a sufficiently large observational sample to test whether random sampling from the IMF is an acceptable 
model in star clusters the available literature was searched for star clusters which are young enough to not have experienced 
supernova events and are dynamically rather un-evolved. For the latter the star clusters should still be e mbedded in their 



natal gas cloud or at least be very young, such that gas expii l sion would not have effected them strongly ijLada et ahl 1984 : 
Goodwin|[l997l : iKroupa et allboOll : [Goodwin fc Bastianlbood : I Weidner et~aLll2007l : IPellerin etallbood : [iBastian fc Goodwin 
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20061 : IWang et aljboosl ). Therefore only clusters younger than 4 Myr have been included in our sample. Additionally, the 



young age limits the amount of mass- loss experienced by massive stars due to stell ar ev olution. 

An allegedly suitable sample of obje c ts dis cussed by Parker fc Goodwin (2007) and Maschberger fc Clarke ( 20081 ) is the 
one compiled by Testi et al. ( 1997 . 19981 . 19991 ) as these authors were explicitly searching for clusters around young A and 
B stars. We do not use the majority of the clusters from these studies for the following reasons: a) the majority are too old 
(> 4 Myr for 25 of 35 objects) or they are b) gas-free. The age limit imposed here is given by the short life time of massive 
stars and to limit stellar mass-loss of the massive stars. Completely gas-free objects are unsuited for the task of this work as 
gas-expulsion will remove large amounts of stars and therefore reduce the mass of the cluster, Med, significantly IjKroupa et al 



200 ll : I Weidner et al.ll2007l ). Th e exception are four obj ects which this sample has in common with the near infra-red study of 



young star- forming reg i ons b y IWang fc Loonevi ( 2007 ) and which are a l l inclu ded in our study. Based on similar arguments 
Maschberger fc Clarkel (|2008l ) also excluded the lTesti etahl (Il997) . |l998| . ll999h sample from their final statistical analysis. 



A very recent additional sample is provided bv lFaustini et al.l (|2009l ). The authors study 26 high-luminosity IRAS sources 
and find that 22 of them show evidence for clustering. They model 9 of these clusters in order to derive cluster masses and 
the mass of the most massive stars. This sample is includ ed in our study , too. But because the results are based on modelling, 
different symbols for them are used in subsequent plots. Faustini et al. ( 20091 ) conclude that the masses of the most-massive 
star in these clusters are also not reconcilable with random sampling of the stars from the IMF. 



2.2 Mass of a cluster versus number of stars in a cluster 



The claim has been made ( Parker fc Goodwin 2007 : Maschberger fc Clarke 20081 ) that the number of stars within a star 
cluster, A'eci, gives a better statistical description of the cluster compared with the cluster mass, Meci, because A^oci is an 
observed quantity and statistically more easily manageable. This is, however, not entirely true as observational biases handicap 
Ncci to a larger extent than Meci. As the lower mass limit of the observations depend on telescope time, distance of the object, 
reddening and observed colour range, the different clusters have to be normalised to the same lower mass limit in order to 
make them comparable. This is done for A'eci in the same way as for Mgci - by extrapolating the stars in the observed mass 
range to a general mass range (0.01 to 150 Mq in this study) with the use of an IMF. Therefore, A'eci is not an observed 
quantity but an estimated one. But the sources for potential error are much larger in the case of Aed than compared with Med, 
as the observed number of stars gives every star the same statistical weight, regardless if it is an M dwarf or an O supergiant. 
But l ow-mass stars and brown dwarfs are easy t o miss due to being faint but also due to un-resolved binarity and crowding of 
stars (jMafz Apellanizll2008l : Iweidner et al.ll2009l ). Very young low-mass PMS stars and brown dwarfs are still difficult to model 
because they are dominated by the unk nown accretion history, and magnetic fields and fast rotation have a strong influence 



( Chabrier et al.l 120071 : iRibas et al.l 120081 ). Therefore is the observational lower mass limit highly model-dependent and has 
large errors. Because the IMF is dominated in number by low-mass objects (85% of all stars are below 0.5 Mq for the IMF 
described in Appendix |X| uncertainties in the lower mass li mit severely affect the Aed estimate. The mass, Med, in contrast 



is far easier to estimate by the number of high-mass stars ( Mafz Apellaniz 20091 ). Likewise, the stellar evolution models of 
massive stars still include large uncertainties. As in the case of low-mass stars, the effects of fast rotation and magnetic fields 
in these stars are not well understood. Massive stars are small in numbers (6% of all stars are above 1 Mq) but dominate the 
cluster in mass (50.7% of the total mass is in stars above 1 Mq for a cluster comprised of 0.01 to 150 Mq objects according 
to the canonical IMF as described in Appendix Because of the intrinsic brightness of these objects they are easy to access 
observationally and difficul t to miss. While the binary frequency might be lower for low-mass stars (~ 35%) compa r ed to 



200' 



i; lApai et al.ll2007l : ISana et al.ll2008l : 



Turner et al. 



2008 



mass i ve stars (20% to 80%.lGarmanv et al.|[l980l: ICarcfa fc MermilliodlboOll: be Becker et ahlbood : iKiminki et aLll2007l : ILucv 



Weidner et al.ll2009l )~ 

for the mass estimate than for the number estimate. If all stars were in un-resolved binaries, Aeci would miss 50% of the stars 



while Med would miss only 16 to 30%, depending on the mass-ratio distribution (jWeidner et al 
to study nimax in dependence of Med rather than Acd. 



20091). We therefore choose 



2.3 Additional issues 



If gas-expulsion already starts early on, before the explosion of supernovae, even the young objects presented here might be 
effected by mass-loss due to the unbinding of stars from the cluster. 

One possible additional effect which might deplete very young star cluster s especially from massive stars are dynamical 

eject ions after stellar encounters in the dense decoupled cores of massive clusters ([Clarke fc Pringlelll992l : iPflamm-Altenburg fc Kroupal 
2006[ ). Unfortunately, this effect is impossible to avoid or to correct for reliably and might lead to an additional underestimation 



of the cluster masses, but it is unlikely th at the most-massive star is ejected from the cluster. 

Already in Weidner fc Kroupal ll2006|) a first se t of young star clusters and their most massive stars were presented. In the 



current contribution the 



Weidner fc Kroupa ( 20061 ) list is included, corrected for a few errors and significantly expanded. The 



sample of 100 both new and previously published star clusters is shown as Tab. IBll in AppendixlB] The table shows two mass 
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Table 1. For these massive stars in star clusters dynamical mass estimates exist from the orbits of binaries. 



Cluster 



Star 



Sp Type 



ynamical 



niold 



' "'IICW 



'■ini new 



Ref. 



Trumpler 14/16 F015B 09.5V 16.0 ± 1.0 23.3 ± 2.0 16.5 ± 1.5 17.9 -3.9/+3.1 (1) 

Trumpler 14/16 F015A 05.5V 30.0 ± 1.0 50.4 ± 6.0 34.2 ± 3.0 37.7 -3.7/+7.3 (1) 

M42 e Orionis CI 06Vpe 35.8 ± 7.2 45.0 ± 5.0 31.7 ± 6.0 34.3 -4.3/+4.7 (2) 

Trumpler 14/16 HD93205A 03V 56.0 ± 4.0 87.6 ± 12.0 58.3 ± 10.0 64.6 -4.6/+5.4 (3) 

Westerlund 2 WR20a A WN6ha 82.7 ±5.5 - - 121.0 -41.0/+29.0 (4) 

NGC3603 NGC 360 3-A1 WN6ha 1 16.0 ±31 120.0 ± 15.0 - 12 1.0 -41.0/+29. (5) 

l lNiemela et al.l l|2006l ). ^ iKraus et aP l|2009f ). ^ iMorrell et~all l|200lh . "^ iNaze et aP l|2008h . ^ ISchnurr etaP l|2008l ) 



values for the mass of the most-massive star. The one in column # 3 is based on the Vacca et al. ( 19961 ) spectral class to stellar 



mass conv ersion. In column # 4, additionally, a new spectral class to mass conversion is used. It is based on lMartins et al.i 
1 20051) and lMartins fc Plej ( 20061 ) who provide two new transformations of 0-star spectral types into masses, which are both 



rather similar. One is based on a theoretical effective temperature scale and the other on an observational one. The authors 
note that their new calibration should represent a significant improvement over previous calibrations, due to the detailed 
treatment of non-LTE line-bl anketinR in their calc ulations. Using the new transformation based on the theoretical effective 
temperature scale (table 1 in Martins et al. 2005h . all the clusters with O star s (m<^16Mp)) in Tab. IBll are re-examined. 
The resulting new spectral masses are corrected for stellar evolutionary effects ( Weidner 20091 ) and the new masses for the 
most-massive stars are compiled in column ^ 4 of the same table. The difference between the old and the new calibration 
is visiualised in Fig. IBll in Appendix [B] As is shown there, in all but f our cases the new ca libration results in stellar masses 
significantly lower than the old values. Note that the new calibration bv lMartins et al.l (120051 1 is only provided up to a spectral 
type of 03. This might not include the most massive stars observed but no general consensus exists in spectral classifying of 
extremely m assive stars. While, t raditionally they would be of spectral type 03 some classify them as spe ctral type 02 o r 
even earlier jWalborn et al.ll2002l ') while others prefere a Wolf-Ra yet star classification ( for example WN 6h, ICrowtherll2007l ). 
In these cases the corrected values based on the [CrowtheiJ (|2007l ) Wolf-Rayet scheme in IWeidnei] (|2009l ) are used. 

An additional complication in the determination of the masses of the most-massive st ars is due to possible binary stellar 
evolution (BSE). Massive stars are often found in close binaries of rather similiar masses llWeidner et al. 200S ) and therefore 



BSE might have a ffected the evolutio n and hence the observational parameters of the stars (jWiiers et al ] |l996l : I Tout et al 
19971 : [Hurlev 20oi : Izhang et allboOSh 



Except for a few cases the cluster masses are derived by extrapolating from the given number of stars above a certain 
mass limit or within certain limits to a mass range of 0.01 to 150.0 Mq with a canonical IMF (see Appendix IX)! . 



Several cluster masses given in lCarpenter et al.l (119931 ) are used as lower limits only in this study because of incomplete- 
ness and uncertain differential reddening. 



In Appendix [C] notes on some individual clusters can be found. 



2.4 Dynamical masses 



In recent years observational techniques allowed to measure masses of very massive stars directly by observing the orbits of 
massive eclipsing binaries. In Table [T] the dynamical mass estimates for six very massive stars are compared with old and new 
spectroscopic estimates. Two of these six stars (WR20a A and O Orionis CI) happen to be the most massive stars in two 
clusters (Westerlund 2 and M42). Also shown in Tab. [1] are the initial masses for these stars, mini new , deriv ed by matching 



the lumi nosity and effe ctive temperature of the newly calibrated O star spectral types by Martins et al. ('2005') with the values 



from the lMevnet fc Ma cder (200S) rotating stellar evolution models for massive stars (for details see I Weid ner 2009,). Generally, 
the new spectroscopic mass estimates from the new calibration agree much better with the dynamical masses than the old 
spectroscopic mass estimates. For the analysis done in this work for WR20a A and O Orionis CI the dynamical masses are 
used for the old calibration and mini new for the new one. 



2.5 The cluster sample 

Table [BT] in Appendex |B] includes the Weidner fc Kroupa ( 20061 ) sample of most- massive stars in star clusters together with 
the new entries compiled here from the literature. Fig. [2] shows the most-massive-star vs. star-cluster-mass relation from this 
table usin g the old stellar ma s ses fo r the O stars. Futhermore, the Fig. shows the theoretical analytic result (the thick solid 
line) from I Weidner fc Kroupal (1200J), which numerically solves eqs.[2]and[3]but with the canonical multi-part power-law IMF 
and assuming a fundamental upper mass limit for stars, mmax* = 150 A4q to arrive at a relation for mmax ~ /ana(Afcci). 
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logio(Meci/Mo) logio(Maci/Me) 
Figure 2. Plot of the most-massive-star vs. star-cluster-mass data using the literature values for the stellar masses and the old effective 
temperature scale calibr ation for 0-stars. Panel A : With error bars and panel B without error bars. The thick solid line represents 
the analytic model from IWeidner fc Kroupal l l2004l . see also Fig. [T] The literature values with ci rcles around them have dynamical 
mass estimates for the most-massive star while the ones with boxes are the sample of lFaustini et all 1^009,3 for which the masses of the 
most-massive stars are only indirectly calculated. The thin solid line marks the identity when a cluster is made of only one star. 



The masses of the most-massive stars derived from the new spectral type to mass conversion are shown in Fig. [Altogether 
with the same lines as in Fig. [5] 



2. 5. 1 Errors 



The error bars for iV/eci in Tab. lBll are either directly taken from the respective literature source, or by assuming an uncertainty 
of 100% in the number of observ ed stars. For mmax the errors are again taken ei ther directly from th e literature or the spectral 
type is converted to mass by the Vacca et al. (199^ tables for column # 3 andlM artins et al" I l|2005l l for column # 4j a-nd the 
spectral subtype -f 1 and -1 is used as the upper and lower limit for the stellar mas^, respectively. The errors in the distance 
and age are from the literature only. 



3 STATISTICAL ANALYSES 

In the supplement of Pflamm-Altenburg fc Kroupa ( 2008h the probability for the ith massive star randomly chosen from a 



number of stars A'^ is given, assuming rrimax* = 150 Mq. For i = 1 (the most-massive star) the probability is, 



N- 



PN(m) yj ^(ni)dmj ^{m), (7) 

with ^(m) (X ^(m) being the probability density distribution and £,(m) the IMF as described in Appendix 1^ 

In order to get the number of stars, A^, required for eq.[7]for a given cluster with Med, an array of cluster masses between 
5 Mq and lO'' Mq is divided by the mean mass, mmcan, of the IMF. For the IMF used here (see appendix |X] for details) 

mmcan = 0.36 Mq, if TUmin = 0.01 Mq and mmax* = 150 Mq. 

Fig. [4] shows the distribution obtained by random sampling (eq. (Tjl for two examples of = 278 {dash-dotted line) and 
A'^ = mil (solid line). For each A' the five statistical values are calculated. 



^ The masses are corrected for stellar evolutionary effects as described in IWeidneil ||2009| ). 

* For example, for an 04 V star the spectral types 05 V and 03 V are used to determine the lower and upper mass limit. 
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logio(Mecl/Mo) loglo( Mecl/Mo) 

Figure 3. Like Fig. [2] but using the new effective temperature-scale calibration for O-stars from [Martins et al. 1 1I2005I') corrected for 
evolutionary effects to initial masses llWeidnedl2009l) . Panel A: With error bars and panel B without error bars. The dynamical mass 
estimates are shown only as open circles as the evolutionary correted initial masses are used for the actual analyses. 
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Figure 4. The solid line shows the distribution of the most-massive star for N = 11111 (Mgd ~ 4000 Mq) and the dash-dotted line is 
the same for = 278 (~ 100 Mq) according to eq. [7] In this plot pjv is normalised to 1 at the most common (mode) value in order to 
give both curves the same height. For the N = 11111 case the following statistical properties are also shown: the vertical dotted line is 
the mode, the vertical dashed line is the median, the vertical solid line is the expectation value and the two vertical long-dashed lines 
mark the l/6th and 5/6th quantiles. See text for further details. 
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• The arithmetic mean or expectation value, marked as a solid vertical line for the N = 11111 case in Fig. [l] is the sum of 
all most-massive stars divided by the number of clusters. 

• The mode value (dotted vertical line in Fig. |4]) marks the most common value (the peak in Fig. |4} of the distribution. 

• The median value (dashed vertical line in Fig.|4]) is the value which divides the distribution in two. 50% of the values are 
above the median while 50% are below. 

• The l/6th quantile (the left long-dashed vertical line in Fig. |4]) is the value below which l/6th of data points lie. 

• The 5/6th quantile (the right long-dashed vertical Ime in Fig. |4| is the value above which 5/6th of data points lie. 

The l/6th and 5/6th quantiles define the region within which lie two thirds of the most- massive stars lie for random 
sampling of stars from the IMF (eq. [7)l . 



3.1 Completeness of the sample 



The completeness of the cluster sample presented here strongly depends on the total number of star clusters expected for 
the Milky Way (MW) whic h are younger tha n 4 Myr. This depends on the assumed current star-formation rate (SFR) of 
the MW (0^ 
1.8 - 2.3 
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yr -. i Diehl et ahlbood and references therein), the slope of the embedded cluster mass function (/9 
Lada fc Ladal 12003 ). where dNcd = M~J^ dMcd is the number of just formed embedded cluste rs with stellar mass 
in th e interval Meci, Mcc\ + dMad, and the assumed lower mass limit for star clusters (5 to 100 Mq, Weidner fc Kroupal 
20061 ') . With the observationally favoured parameters being SFR = 4.0 Mq yr~^, P — 2.0 and M, 



yr ', p ^ -z.u ana iWeci min = 5Mq. The total 
number of young star clusters in the MW lies therefore between 10* and 10® clusters. The majority of these have masses less 
than 100 Mq and any surveys of them are severely incomplete. For a completeness estimate we therefore restrict ourselves to 
clusters more massive than 1000 Mq as they are far fewer in numbers and more easily identified in the MW. For the whole 
range of parameters of the MW the number of young star clusters more massive than 1000 Mq lies somewhere between 160 
and 4452, with 1478 being the value for the observationally favoured parameters. The sample shown in Tab. IBll includes 30 
(-5/-I-6) clusters which are in the MW and more massive than 1000 Mq within the uncertainties. This suggests that between 
18.8% (-3.1/4-3.7) and 0.7% (± 0.1) of all such clusters are in the sample, with 2.0% (-0.3/-I-0.4) for the favoured parameters. 
Therefore, one has to keep in mind that any statistical results are possibly limited by the incompleteness of the cluster sample. 



3.2 Statistical tests 

In panel A of Fig. [5]the mode, mean, median and l/6th and 5/6th quantiles for a fundamental upper mass limit of mmax* ~ 
150 Mq are shown together with the data points from column # 4 and the clusters from column # 3 from Tab. IBll which 
have not been changed by the re-calibration. 

Three different statistical tests are applied to the data in order to verify whether or not the observed most-massive stars 
are consistent with being randomly drawn from the IMF. 



3.2.1 Percentage of stars between the l/6th and 5 /6th quantiles 

As is visible in this Fig. [5] there is a general change in behaviour of the data points around a cluster mass of about 100 Mq 
and around 1000 Mq with respect to what is expected from random sampling. Below the 100 Mq limit the data show a 
larger spread while above 1000 Mq the slope of the rrimax-Afcci-relation changes. Panel A of Fig. [6] shows the percentage of 
the most-massive stars within the l/6th and 5/6th quantiles in three samples, one for the clusters below 100 Mq, one for the 
clusters between 100 and 1000 Mq and one for the ones above 1000 Mq . Additionally, the figure shows the same numbers for 
different assumptions on the fundamental upper mass limit (mmax*) for stars. Here the clusters above 1000 Mq (filled and 
open triangles) are far below the 2/3rd range which would be expected from random sampling. The clusters below 100 Mq 
(filled and open circles) and the intermediate clusters (100 to 1000 Mq, filled and open squares) are very tightly within the 
l/6th and 5/6th quantiles. About 90% and 78% of the clusters are within the range, respectively. In panel B of Fig. |6]the 
same is shown but including the error bars for mmax and Mod from Tab. IBll by making the same calculations as before but 
using the minimal and maximal values for mmax and A/oci- The low-mass clusters are still more tightly distributed within the 
l/6th and the 5/6th quantiles than expected. The intermediate and high- mass cluster seem to be consistent with random 
sampling when the maximum effect of the errors is applied to the data. 



3.2.2 Distribution around the Median 



Also important is the distribution of the mmax obs values around the median of the expected distribution for random sampling. 
The median is the statistical value for which 50% of th e data should l ie abo ve and below. For the whole sample 25.7% are 
above the median and 74.3% below if one uses the new iMartins et al] l|2005ll O star mass scale and assumes a fundamental 
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Figure 5. Panel A: Most-massive star vs cluster mass. The dots are the observed values from column # 4 from Tab. IBll The two open 
ci rcles indicate existing dynamical estimates for the present-day mass of the most-massive stars. The boxed data are from the sample 
of iFaustini et al"] | |2009| ). The dotted line refers to the mode value for random sampling, the short-dashed line to the median value, the 
curved solid line marks the mean value and the two long-dashed lines are the l/6th (lower) and 5/6th (upper) quantiles between which 
2/3rd of the data points should lie if they were randomly sampled from the IMF over the mass range 0.01 Mq to mmax* = 150 Mq. The 
thin solid line to the left marks the identity where the cluster is made-up only of one star. Panel B: The same as panel A but assuming 
a fundamental upper mass limit for stars of mmax* = 50 Mq instead of 150 Mq . 



upper mass limit of mmax* = 150 Mq. In the sub-sample of clusters below 100 Mq there are 56.8% above and 43.2% below the 
median, for the clusters with 100 < Mod ^ IOOOMq there are 13.3% above and 86.7% below the median while for the high-mass 
clusters 2.9% and 97.1% are above and below the median, respectively. In Fig.[7]the distribution of mmax median — mmax obs new 
is shown for the whole cluster sample {panel A) and for the clusters below 10^ Af© (panel B). For the old O star mass scale the 
distribution is shown in Fig.|8l The percentages in the case of the old O star mass scale are 35.6/64.4% for the total sample and 
56.8/43.2%, 40.0/60.0% and 8.8/91.2% for, respectively, the clusters below 100 Mq, the clusters with 100 < A'/ed < IOOOMq 
and the clusters above 10^ Mq. 



3.2.3 Wilcoxon-Signed-Rank-Test 

The Wilcoxon signed rank test ( Bhattacharvva fc Johnson 1977 fl tests whether or not the data is consistent with being 
symmetrically distributed around the median. It reveals for the new calibration a probabilitjjf], p(Mcci ^ 100A/q), of 0.014 for 
clusters with masses smaller or equal to 100 Mq, a p(100Mq < Moei < IOOOMq) of 1.9 • 10"'' for cluster masses between 100 
and 1000 A^q and a p(Mcci > IOOOA/q) of 2.8 • 10"® for the clusters above 1000 Mq. For the old calibrations the probabilities 
are p(A/ed sS IOOA/q) = 0.014, p(100Mq < A/ed < IOOOMq) = 0.035 and p(Mod > IOOOMq) = 1.2 ■ 10"**. 



3.3 Dependence on the high-mass IMF slope 

The general assumption in this paper, that the stars in a star cluster follow a universal IMF which is characterised by a 
Salpeter/Massey slope of 2.35 for all stars above 0.5 Mq, is strongly supported by almost all observational evidence (see 
appendix |X] for a list of references). However, if the IMF-slope for hig h-mass stars is steepe r than 2.35, it is not sure whether 



or not a fundemental upper mass limit exists, as was pointed out by Oev fc Clarke! (2005). Fig. |9] shows the mean, median 



mode, l/6th and 5/6th quantiles for two different assumptions of the high-mass slope of the IMF. In Panel A the slope for 
stars more massive than 25 Mq is changed to Qa = 3.0 and in Panel B to as = 4.1, while Q2 ~ 2.35 for stars between 0.5 

^ A short introduction into the test and pre-calculated tables for the probabilities for different A'^ can be found at: 
|http: // comp9 . psych. Cornell, edu /Darlington / index.htm 

° ~ 0.1 is the highest probability the Wilcoxon Signed Rank test allows for. 
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Figure 6. Percentage of most-massive stars within the l/6th to 5/6th quantiles for three different assumptions of the fundamental upper 
mass hmit, rrimax*, of 50, 100 and 150 Mq. Solid symbols refer to the new O star caUbration in Tab. IBll and open symbols to the old 
one. Circles mark the percentages for clusters below IO^Mq while squares are for 100 < Mgd ^ IO^A/q and triangles mark clusters 
more massive than 10^ Mq. Panel A: These values do not take into account the errors in Mgci and mmax in Tab. IBll while in panel B 
these errors are included. The horizontal dotted line at 66% marks the fraction of stars that ought to lie between the l/6th and 5/6th 
quantiles if random sampling between 0.01 Mq and rrimax* were true. 




3 4 

logio(Meci/Mo) logio(Mee,/Me) 
Fig ure 7. Distanc e of th e observed most-massive star from the expected median for random sampling up to rrtmax* = 150 Mq. Here the 
new lMartins et al. I l l2005D mass scale is used for O stars. Panel A show the whole sample while panel B contains only the clusters below 
IO'^Mq. For random sampling over the mass range 0.01 Mq to 150 Mq the data ought to be distributed symetrically about Ammax = 0. 
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and 25 A4q in both cases. 25 Mq was chosen because it is the mmax value of the plateau like feature for clusters with masses 
between 10"^ and 4 • 10^ Mq. Only for a slope as steep as as = 4.1 are 50% of stars above and below the median. A slope 
steeper than as = 3.7 is needed in order to have more than 60% of the stars within the l/6th and 5/6th quantil es. Such steep 
slopes for the h i gh-mass IMF within star clu sters are clearly ruled out by the current state of observations ( Massev 19981 : 
Elmegreen|[l999l : lKroupall200ll : lLarsonll2002bl ). 
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Figure 10. Like Fig. |3]but only tiie clusters wliicli are in common with the lPfalzneJ ||2009| ) sample. Filled circles are clusters with radii 
smaller than 1 pc while the open circles are the clusters with larger radii. 



3.4 Dependence on the environment 

Very recently, Pfalznei ( 20091 ) studied the dissolution behaviour of young (1 to 20 Myr) massive star clusters (2000 to 50000 



Mq). She found that her sample of 23 clusters can be divided into two groups, loose clusters {Red > 1 pc) and tight clusters 
(J?eci < 1 pc), where i?oci is her estimated cluster radius. The radii of the groups each follow a rather tight sequence with 
time. While the tight clusters expand from ~ 0.5 pc to 3 pc the loose ones evolve from 4 pc to 20 pc on the same time scale, 
parallel to the tight ones. Of these 23 clusters 10 are included in our cluster sample. 5 of them are tight clusters ([OBS 2003] 
179, Westerlund 2, NGC 3603, Trumpler 14, Arches) and 5 are loose ones (NGC 7380, NGC 2244, IC 1805, NGC 6611, Cyg 
0B2) . When comparing the most- massive stars against the cluster mass of these two subsets, as is done in Fig. 1101 it seems 
that the clusters which form tighter and are therefore more dense, have on average a more massive maximal star, while the 
loose clusters prefer less massive maximal stars. For the tight subset a linear function can be fitted with a slope of 0.09 ± 0.39 
and a rather low linear correlation coefficient of 0.35 ± 0.47. The slope for the loose sample is 0.27 ± 0.16, somewhat steeper 
than for the tight sample within the error bars, but the linear correlation coefficient is much larger, about 0.92 ± 0.08 . The 
combined sample has a slope of 0.22 ± 0.23 with a linear correlation coefficient of 0.52 ± 0.47. Therefore, the difi'erence in 
slopes might be indicating a physical dependence of the mass of the most massive star not only on the cluster mass (previous 
sections) but also on the cluster density. But the large error bars make a more definite statement difficult. Also it should 
be noted here that the R^ci estimates for the all the loose clusters of the Pfalzner 1 2009l )-sample are the measured median 
distances of early B type stars in theses clusters (iWoffl et a l. 2007) and therefore might not be directly comparable to radii 

arrived at with difi'erent methods. 

Low-mass young clusters arc found to be generally small ( 1 pc, Testi et al. 19981 : Gutermuth et al. 20051 : Testor et al 



[2005; Rathbornc et al. 2006; Schccpmakor ct al. 2007) so it is unclear if and how such a correlation between the most-massive 
star and the cluster radius extends to lower masses. 



3.5 A simple Model 

A simple yet sufficient model to describe the plateau of most-massive stars between 1000 and 4000 Mq and the behaviour at 
higher cluster mass might be the following. The model assumes that the mass of the most-massive star is linked to the proto 
cluster mass due to stellar feedback. For a range of cluster masses, Med (10 to 10® Mq), cluster radii, i?cci (from 0.1 to 1.0 
pc), and star-formation efficiencies, SFE (0.3 to 0.8), the velocity dispersion, a, is calculated by 



2 iiociSFE' 



(8) 
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Figure 11. The dashed lines show the dependence of the velocity dispersion, a, in the cluster on Mg^i for a range of cluster radii (0.1 
to 1 pc) and star-formation efficiencies (0.3 to 0.8). The horizontal lines mark the typical range of velocities for ionised gas, V[an (10 to 
20 kms~^), while the vertical lines are the plateau of most-massive stars for M^ci from 1000 to 4000 Mq. In the shaded region are lying 
models for which a is larger than Hion and Med is larger than 1000 Mq. 



where G is Newton's gravitational constant and SFE = — -, with Mgas being the residual gas mass in the cluster 
forming volume. 

Fig. [11] shows a within the proto cluster as a function of Med- It is compared with the typical velocity of ionised gas, 
Uion, which is about 10 to 20 kms~^. As is visible in Fig. 1111 a is larger than «ion for clusters with masses larger than a couple 
of hundred Mq , regardless of the radii and SFEs. Therefore, it seems possible that such clusters are able to retain the ionised 
gas longer - allowing the stars to accrete further mass. The fact that a already overcomes Hion at rather low Med for small 
Red and low SFE can be seen as an indication that low-mass clusters might have lower SFEs than massive clusters. 



4 RESULTS & DISCUSSION 

We have studied the possible dependence of the mass of the most-massive star, mmax, on the stellar mass, Med, of the host 
birth cluster. To this effect we have significantly increased the data sample mmax(M ed). 



Using the new spectral-type-stellar-mass conversion from iMartins et al.l (12005! ) and the here presented sample of most 



massive stars in star clusters, it has been shown here that the observed sample divides into three sub-samples, the first being 
clusters with M^d < 100 A'Iq, followed by clusters between 100 and 1000 Mq and clusters with Med > 1000 Mq. Furthermore, 
there is a plateau of constant rrimax ~ 25 Mq for clusters with masses between 1000 and 4000 Mq. 

• A'fed < IOOMq: The percentage of stars between the l/6th and 5/6th quantiles is 89% (83% when taking the error bars 
into account) which is too tight for random sampling (66%). Such a distribution is highly unlikely with a chance of only 0.2% 
(0.1% with errors) when calculated from a Binomial distribution. But the distribution around the median and the Wilcoxon 
singed rank test are compatible with random sampling at a significance of 2 percent. 

• 100 < Med ^ IOOOMq: 77% (70% with errors) of the stars are within the l/6th and 5/6th quantiles which is somewhat 
tighter than expected for random sampling (66%). The probability of this to occur is rather high with 8% (13% with errors). 
But 87% of all clusters are below the random- sampling median where only 50% would be expected and the Wilcoxon singed 
rank test gives a very low probability (1.9 • 10~^) that the data are distributed symmetrically around the median. 

• Med > IOOOM0: Only 12% of the data points (66% with errors) are in the 2/3rd interval which is far below the 
expectation from random sampling (66%). The probability for a random occurance of such a low number with the 2/3rd 
interval is 4 ■ 10~^^. Furthermore 97% of the data points are lower than the median and the Wilcoxon singed rank test results 
in a very low probability (2.8 ■ 10~^) for a symmetric distribution, too. 



The clusters in the mass range below 100 Mq are the ones most compatible with the hypothesis of being randomly 
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Figur e 12. The mass of the most-massive star vs the mean mass of the so-called "superclusters" constructed as in lSelman fc Melnickl 
KOOSh . It shows whether or not the observed clusters of higher masses can be made by adding up large numbers of low- mass clusters. 
The obvious trend of increasing rrtmax with < M^ci > is a clear indication against such a conclusion. 



sampled from the IMF. This is also roughly the range of clusters studied by Maschberger fc Clarkel ( 20081 ) . Their result, that 
the most-massive stars in these clusters could be randomly drawn from a universal IMF, is therefore in accordance with 
our conclusions. The difference is that here it is shown that this assumption can not be generalised for more massive/richer 
clus ters. 

Selman fc Melnick 1 20081 ) argue that the claim reached by Weidner fc Kroupal 1I2OO6I ). that there exists a minax(Afeci)- 



relation, is due to a size-of-sample effect in the data used by I Weidner fc Kroupal (|2006l ). We now apply their analyses to our 
new data set. In appendix A of their paper they use a method of adding up some clusters of the Weidner fc Kroupa ( 20061 ) 
sample to so-called "superclusters" of the same mass as NGC 6530 (about 1000 Mq in the new sample presented here). By 
comparing the mean mass of the synthetic superclusters with the most-massive star of the component clusters, they show 
that there is no trend for the most-massive stars to be more massive with cluster mass. Here we repeat the same method 
with our new sample of clusters. AH possible combinations to reach the mass of NGC 6530 from within the sample are used 
and the most-massive star is plotted over the mean cluster mass, < M e d >, in Fig. 1121 As is see n in the figure the mass o f 
the most-massive star increases with < Med >• The Selman fc Melnick ( 20081 ) explanation for the Weidner fc Kroupa ( 20061 ) 
result therefore fails for the new sample. 

These results strongly suggest an underlying physical mmax-Afcd-relation. They contradict the hypothesis that star 
clusters are populated with stars by random sampling from the IMF. Only when taking into account the full range of the 
error bars and a very unlikely low fundamental upper mass limit of mmax* = 50 Mq would the complete sample mostly agree 
with random sampling. But in such a case no stars above 50 Mq would exist, a result clearly disproved by the dynamical 
mass measurements for the massive stars in Westerlund 2 and NGC 3603 (see Tab.[T}. 

The general trend of the most-massive star with cluster mass and the observed plateau between the two cluster mass 
regimes is therefore most likely a general result of the star-formation process within cluster-forming molecular cloud cores. 
Several different mechanisms might be responsible for the non-random behaviour of the formation of the most-massive star 
in star clusters. One such model is explored in § 13.51 where the velocity dispersion within the cluster-forming cloud core is 
used as a measure for the binding energy of the cloud, and is compared with typical velocities of ionised gas which acts as a 
proxy for the radiative feedback of the stars. This simple model is already in qualitative agreement with the data, but more 
detailed studies of how the radiative and mechanical feedback of massive stars scales differently than the binding energy are 
needed. This may result in a critical Afed limit at which the one dominates over the other. 

Another possible explanation for the existence of an mmax (Med )-relation might be given by dry mergers. In this scenario 
massive stars form in smaller sub-clusters which are quickly evacuated by their feedback and these sub-clusters then merge 
nearly gas-free, allowing only for very little additional accretion, ie mass growth. Only for initially very massive giant molecular 
clouds, more gas might be accreted during and after the merging of the sub-clusters. 
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The interesting split of the massive clusters into a tight and a loose subset by IPfalznerl (|2009l . see § I3.4|l can be used 
as an additional constraint on the mmax-Afod-relation. The loose cluster stars form predominately by free-fall collapse of 
dense cores with little or no further gas accretion into the cluster. But for the tight (high-density) clusters cluster-potential- 
assisted-accretion is possible which allows for more massive stars to form in these objects. Also stellar collisions, mergers and 
competitive accretion might play a role in these dense clusters. 

A more detailed study of the possible mechanisms to explain the here presented observational evidence for a physical 
relation between mmax and Meci will be presented in a follow-on paper. 

As the high-mass regime is most important for the question whether the integrated IMF of a galaxy is similar to the IMF 
derived locally on star cluster scales or not, this discardation of random sampling naturally leads the IGIMF being steeper 
than expected from individual star clusters. Since the majority of stars seem to form in star clusters but also these clusters are 
distributed according to a mass function which is dominated by lower mass clusters, the apparent non-randomness of these 
clusters lead to fewer OB stars per star in a galaxy than expected from random samplinj^l- 
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APPENDIX A: THE STELLAR INITIAL MASS FUNCTION 



The following multi-component power-law IMF is used throughout the paper: 



^(m) = k < 



\— ) ,mu < mo, 

[ — ) ,mo<m<mi, 



(Al) 



with exponents 

ao = +0.30 
ai = +1.30 
«2 = +2.35 
03 = +2.35 



0.01 ^ m/M0 < 0.08, 
0.08 ^ m/M0 < 0.50, 
0.50 ^ m/M© < 1.00, 
1.00 < m/MQ. < mma 



(A2) 



where dN = ^{m)dm is the number of stars in the mass interval m to m + dm. The exponents Qi represent the standard 
or canonical IMF an d have been corrected for unresolved multiple systems I Kroupa 2001 . 20021 : Thies fc Kroupa 2007 . 20081 : 
Weidner et al ] |2009l ). The advantage of such a multi-part power-law description is the easy integrability and, more importantly, 
that different parts of the IMF can be changed readily without affecting other parts. Note that this form is a two-part power-law 
in the stellar regime, and that brown dwarfs contribute about 4 per cent by mass only and need to be treated as a separate 
popu lation such that the IMF has a discontinuity near mn = 0.08 M© with A;' ~ i jXroupa et al.ll2003l:lThies fc KroupalboO?! , 
20081 ) . A log- normal form below 1 M© wi th a power-law extension t o high masses was suggested by Chabriei (2003) but is 
indistinguishable from the canonical form (jDabringhausen et al ]l2008') and does not cater for the discontinuity . The canonical 
IMF is today understood to be an invariant Salpeter/Massey power-law slope ( Salpeteil 



bein^ ^ indep endent of the cluster densi ty ari d met allicity for metallicities 0.002 



ty^ 

2000, 2002; ^ Parker et allboOll: 'Masscv''l998l. |2002| . I2OO3I : iLarsonl [iooial lbl: IWvse et al 
2004 : .Pflamm-Ahenburg fc Kroupa, 2009 1 



1955 



Massev 



^assey fc Hunterlll998 



2003) above 0.5 M, 



2002 



BeU et al.ll2003 



Sirianni et al 



Piskunov et al 



The basic assumption underlying our approach is the notion tha t all stars in every cluster are drawn from this same 
universal parent IMF, which is consistent with observational evidence ( Elmegreen 19991 : Kroupa 2001). 

It should be n oted here that, while not indicated in eq. IA2I there is evidence of a maxim al mass for stars (rrtmax ^ 
?Ttmax* ~ 150 M © , Weidner fc Kroupa 2004 ). a result confirmed by several independent studies ( Oev fc Clarkel 2005 : Figei 
2005l : lKoen|[2006l ). 



^ The IGIMF concept is discussed in great detail in iKroupa fc Weidneil ||2003| ) and IWeidner fc Kroupal i2Q0^ and subsequent papers. 
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APPENDIX B: THE CLUSTER SAMPLE 



Table Bl. Literature data for the empirical cluster masses (Med), maximal star masses (mmax.obs) within these clusters, cluster ages (age), distances (D), the numbers of stars above 
or within certain mass limits (in Mq), the name and the spectral type of the most-massive star and references for the data. 



Designation 


Mecl 


^max obs old 


'^max obs new 


age 


D 


# of stars 


Id mmax 


Sp Type 


Ref. 




[Mq] 


[Mq\ 




[Myr] 


[pc] 










IRAS 05274+3345 


14 -7/+15 


7.0 ± 2.5 




1.0 


1800 


15 > 0.24 


- 


B2 


(1) 


Mol 139 


16 + 8 


2.9 ± 2.0 




<1 


7300 


- 


- 


- 


(2) 


Mol 143 


21 ± 10 


3.1 ± 2.0 




<1 


5000 


- 


- 


- 


(2) 


IRAS 06308+0402 


24 -13/+25 


11.0 ± 4.0 




1.0 


1600 


16 > 0.37 


- 


B0.5 


(1) 


VV Ser 


25 -13/+27 


3.3 ± 1.0 




0.6 


440 


24 > 0.3 


VV Ser 


B9e 


(3) 


VY Mon 


28 -15/+29 


4.1 ± 1.0 




0.1 


800 


26 > 0.3 


VY Mon 


B8e 


(3) 


Mol 8A 


30 ± 15 


3.8 ± 2.0 




<1 


11500 


- 


- 


- 


(2) 


IRAS 05377+3548 


30 -15/+32 


9.5 ± 2.5 




1.0 


1800 


31 > 0.24 


- 


Bl 


(1) 


Ser SVS2 


31 -16/+31 


2.2 ± 0.2 




2.0 


259 ± 37 


50 > 0.17 


BD 01° 3689 


AO 


(4) 


Tau-Aur 


31 -18/+36 


7.0 ± 2.5 




1-2 


140 


123 > 0.02 


HK Tau/Gl 


B2 


(5) 


IRAS 05553+1631 


31 -16/+33 


9.5 ± 2.5 




1.0 


2000 


28 > 0.28 


- 


Bl 


(1) 


IRAS 05490+2658 


33 -17/+36 


7.0 ± 2.5 




1.0 


2100 


30 > 0.29 


- 


B2 


(1) 


IRAS 03064+5638 


33 -17/+36 


11.0 ± 4.0 




1.0 


2200 


27 > 0.31 


- 


B0.5 


(1) 


IRAS 06155+2319 


34 -18/+35 


9.5 ± 2.5 




1.0 


1600 


38 > 0.21 


- 


Bl 


(1) 


Mol 50 


36 ± 18 


3.5 ± 2.0 




<1 


4900 


- 


- 


- 


(2) 


Mol 11 


47 ± 20 


3.8 ± 2.0 




<1 


2100 


- 


- 


- 


(2) 


IRAS 06058+2138 


51 -27/+54 


7.0 ± 2.5 




1.0 


2000 


26 > 0.49 


- 


B2 


(1) 


NGC 2023 


55 -28/+58 


8.0 ± 2.0 




3.0 


400 


21 > 0.6 


HD 37903 


B1.5V 


(6) 


Mol 3 


61 ± 20 


3.7 ± 2.0 




<1 


2170 


- 


- 


- 


(2) 


Mol 160 


63 ± 20 


4.3 ± 2.0 




<1 


5000 


- 


- 


- 


(2) 


NGC 7129 


63 -33/+104 


9.2 ± 3.0 




0.1 


1000 


53 > 0.3 / 3 > 3 


BD 65° 1637 


B3c 


(7) 


IRAS 06068+2030 


67 -35/+70 


11.0 ± 4.0 




1.0 


2000 


59 > 0.28 


- 


B0.5 


(1) 


IRAS 00494+5617 


71 -37/+74 


9.5 + 2.5 




1.0 


2200 


58 > 0.31 


- 


Bl 


(1) 


IRAS 05197+3355 


72 -38/+75 


11.0 ± 4.0 




1.0 


3200 


34 > 0.50 


- 


B0.5 


(1) 


Cha I 


80 -46/+91 


5.0 ± 3.0 




2.0 


170 ± 10 


237 > 0.04 


HD 96675 


B6IV/V 


(8) 


V921 Sco 


71 -36/+429 


14.0 ± 4.0 




0.1-1 


800 


33 > 0.5 


V921 Sco 


BOe 


(7) 


IRAS 05375+3540 


73 -38/+78 


11.0 ± 4.0 




1.0 


1800 


74 > 0.24 


- 


B0.5 


(1) 


IRAS 02593+6016 


78 -41/+81 


15.0 ± 5.0 




1.0 


2200 


61 > 0.31 


- 


BO 


(1) 


Mol 103 


80 ± 20 


4.0 ± 2.0 




<1 


4100 


- 


- 


- 


(2) 


NGC 2071 


80 -44/+89 


4.0 ± 2.0 




1.0 


400 


105 > 0.2 


V1380 Ori 


B5 


(9) 


MWC 297 


85 ±60 


8.3 -1.3/+13.7 


8.3 -1.3/+13.7 


0.1-1 


250 or 450 


24 > 0.3 


- 


B1.5V or 09c 


(3) 


IC 348 


89 -46/+92 


6.0 ± 1.0 




1.3 


310 


173 > 0.1 


BD 31°643 


B5V 


(10) 


BD 40° 4124 


90 -49/+106 


12.9 -6.0/+ 2.0 




0.1-6 


1000 


74 > 0.3 / 3 > 3 


BD 40° 4124 


B2c 


(3) 


p Oph 


91 -46/+93 


8.0 ± 1.0 




0.1-1 


130 ± 20 


78 > 0.3 


p Oph 


BIV 


(11) 


IRAS 06056+2131 


92 -49/+97 


7.0 ± 2.5 




1.0 


2000 


85 > 0.28 


- 


B2 


(1) 


IRAS 05100+3723 


98 -51/+103 


15.0 ± 5.0 




1.0 


2600 


63 > 0.38 


- 


BO 


(1) 


R CrA 


105 -55/+114 


4.0 ± 2.0 




1.0 


130 


55 > 0.5 


R CrA 


A5eII 


(12) 


NGC 1333 


105 -54/+111 


5.0 ± 1.0 




1-3 


250 


134 > 0.2 


SSV 13 




(13) 


Mol 28 


105 ± 20 


9.9 ± 2.0 




<1 


4500 








(2) 


IRAS 02575+6017 


111 -57/+116 


9.5 ± 2.5 




1.0 


2200 


91 > 0.31 




Bl 


(1) 


W40 


144 -80/+576 100 


10.0 ± 5.0 




1-2 


600 


3 > 4 


IRS 2a 




(14) 


a Ori 


150 -76/+155 


20.0 ± 4.0 


17.3 -4.3/+4.7 


2.5 


360 ± 60 


140 ± 10 (0.2-1.0) 


a- Ori A 


09-9.5V 


(15) 


NGC 2068 


151 -86/+169 


5.0 ± 3.0 




1.0 


400 


192 > 0.2 


HD 38563A 


B4V 


(16) 


NGC 2384 


189 -95/+192 


16.5 ± 1.5 




1.0 


2100 


7 > 3 


HD 58509 


B0.5III 


(17) 



Designation Afcd "Jmax obs old "'max obs new ^§'2 ^ of stars Id Tflraax Sp Type Ref. 
[Mq\ [Mq] [Mq] ^vlyij [pc] 



Mon R2 


225 -117/+236 


15.0 ± 5.0 




0-3 


830 ± 50 


309 > 0.15 


IRS ISW 


BO 


(18) 


IRAS 06073+1249 


239 -120/+242 


11.0 ± 4.0 




1.0 


4800 


25 > 1.47 


- 


B0.5 


(1) 


Trumpler 24 


251 -131/+291 


14.5 ± 2.5 




1.0 


1140 


4 > 5 


GSC 7872-1609 


WN 


(19) 


IC 5146 


293 -226/+305 


14.0 ± 4.0 




1.0 


900 


238 > 0.3 / 5 > 3 


BD 46° 3474 


BOe 


(20) 


HD 52266 


400 ± 350 


28.0 ± 3.5 


20.0 -5.0/+7.0 


<3.0 


1700 ± 1000 


4 ± 2 > 4 


HD 52266 


08-9V 


(21) 


HD 57682 


400 ± 350 


28.0 ± 3.5 


20.0 -5.0/+7.0 


<3.0 


? 


4 ± 5 > 4 


HD 57682 


08-9V 


(21) 


Alicante 5 


461 -234/+516 


12.0 ± 4.0 




<3.0 


3600 +600/-400 


22 > 2.5 


A47 


B0.7V 


(22) 


Ccp OB3b 


485 -243/+497 


37.7 ± 5.0 


27.2 -5.2/+8.8 


3.0 


800 ± 100 


12 > 4 


HD 217086 


07Vn 


(23) 


HD 153426 


500 ± 350 


40.0 ± 6.5 


28.9 -6.9/+7.1 


<3.0 


7 


5 ± 4 > 4 


HD 153426 


06.5-7V 


(21) 


NGC 2264 


525 -267/+537 


25.0 ± 5.0 


27.2 -5.2/+8.8 


3.0 


760 ± 100 


1000 > 0.08 


S Mon / HD 47839 


07Ve 


(24) 


Sh2-294 


525 -267/+540 


12.5 ± 2.5 




4.0 


3200 


155 > 0.7 


S294B0.5V 


B0.5V 


(25) 


RCW 116B 


536 -276/557 


21.0 ± 5.0 


16.9 -5.9/+7.1 


2.5 


1100 


102 > 0.95 


- 


B1V-08V 


(26) 


NGC 6383 


561 -281/+563 


37.7 ± 5.0 


27.2 -5.2/+8.8 


2.0 


1300 ± 100 


21 > 3 


HD 159176 


07V -t- 07V 


(27) 


Alicante 1 


577 -290/+583 


45.0 ± 5.0 


34.3 -8.3/+10.7 


2-3 


4000 ± 400 


38 > 2.0 


BD 56° 864 


06V 


(28) 


HD 52533 


621 -417/+1077 


26.7 ± 3.0 


19.1 -4.1/+2.9 


<3.0 


? 


15 ± 5 > 4 


HD 52533 


08.5-9V 


(21) 


Sh2-128 


666 -342/+736 


37.7 ± 5.0 


27.2 -5.2/+8.8 


2.0 


9400 


7 > 7 


- 


07V 


(29) 


NGC 2024 


690 -350/+706 


20.0 ± 4.0 


15.3 -6.8/+8.7 


0.5 


400 


309 > 0.5 


IRS 2b 


08V - B2V 


(30) 


HD 195592 


725 -364/+757 


40.0 ±10.0 


32.4 -6.4/+6.6 


<3.0 


7 


18 ± 3 > 4 


HD 195592 


06-6.5V 


(21) 


Sh2-173 


748 -395/+901 


25.4 ± 5.0 


18.3 -4.3/+3.7 


0.6-1.0 


2500 ± 500 


7 > 7 


BD 60° 39 


09V 


(31) 


DBSB 48 


792 -416/+1126 


56.6 ± 15.0 


42.3 -8.3/+12.7 


1.1 


5000 ± 700 


5 > 10 


- 


05V 


(32) 


NGC 2362 


809 -409/+823 


43.0 ± 7.0 


35.3 -5.3/+34.7 


3.0 


1390 ± 200 


353 > 0.5 


T CMa 


09Ib 


(33) 


Pismis 11 


896 -448/+938 


40.0 -0.0/+40.0 




3-5 


3600 -t-600/-400 


43 > 2.5 


HD 80077 


B2Ia 


(22) 


NGC 6530 


1075 -545/+1097 


56.6 ± 11.0 


30.5 -6.5/+8.5 


2.3 


1350 ± 200 


620 > 0.4 


HD 165052 


06.5V 


(34) 


FSR 1530 


1410 -707/+1581 


30.0 ± 15.0 




<4.0 


2750 ± 750 


35 > 4 


[M81]I-296 


- 


(35) 


Berkeley 86 


1440 -730/+1470 


40.0 ± 8.0 


21.7 -4.7/+7.3 


3-4 


1700 


340 > 0.8 


HD 193595 


08V(f) 


(36) 


NGC 637 


1682 -854/+1726 


30.8 ± 15.0 


21.7 -6.7/+11.3 


4.0 


2160 


583 Mq > 1.6 


- 


^ 08 


(37) 


W5Wb 


1734 -874/+1757 


34.1 ± 5.0 


24.1 -4.1/+8.9 


2.0 


2000 


300 > 1 


BD 60° 586 


07.5V 


(38) 


Stoek 16 


1857 -955/+2045 


43.0 ± 10.0 


32.1 -7.1/+7.9 


4.0 


1650 


16 > 8 


HD 115454 


07.5ni 


(39) 


ONC 


2124 -1078/+2175 


35.8" ± 7.2 


34.3 - 5.7/+10.7 


<1.0 


414 ± 7 


3500 (0.1 - 30) 


© Orionis CI 


06Vpe 


(40) 


RCW 38 


2251 -1132/+2276 


50.4 ± 10.0 


37.7 -7.7/+12.3 


<1.0 


1700 


2000 > 0.25 


IRS 2 


05.5V 


(41) 


Bochum 2 


2284 -1302/+3523 


37.7 ± 4.0 


27.2 -5.2/+8.8 


2-4 


2700 


4 > 16 


BD 00° 1617B 


07V 


(42) 


Berkeley 59 


2310 -1168/+2417 


37.7 ± 4.0 


27.2 -5.2/+8.8 


2.0 


1000 


41 > 5 


BD 66° 1675 


07V 


(43) 


IC 1590 


2376 -1245/+2799 


41.0 ± 5.2 


30.5 -6.5/+8.5 


3.5 


2900 


14 > 10 


BD 55° 191 


06.5V 


(44) 


W5E 


2614 -1323/+2667 


37.7 ± 5.0 


27.2 -5.2/+8.8 


2.0 


2000 


400 - 500 > 1 


HD 18326 


07V 


(38) 


W5Wa 


2651 -1338/+2690 


52.0 ± 10.0 


37.3 -6.3/+7.7 


2.0 


2000 


400 - 500 > 1 


HD 17505 


06.5HI 


(38) 


NGC 1931 


3128 -1564/+3163 


30.8 ± 15.0 


21.7 -6.7/+11.4 


4.0 


3086 


848 Mq > 2.39 


- 


RJ 08 


(37) 


LH 118 


3746 -1918/+4077 


56.6 ± 7.0 


42.3 -8.3/+12.7 


3.0 


48500 


28 > 9 


LH 118-241 


05V 


(45) 


NGC 2103 


3853 -1937/+3905 


80.0 ± 20.0 


89.3 -24.3/+30.7 


1.0 


48500 


26 > 10 


Sk -71° 51 


02V({P)) 


(46) 


NGC 7380 


4527 -2290/+4611 


47.8 ± 5.0 


36.0 -6.0/+9.0 


2.0 


3700 


42 (6 - 12) 


HD 215835 


05.5-6V((f)) 


(47) 


NGC 6231 


4595 -2312/+4676 


59.0 ± 10.0 


51.9 -20.9/+18.1 


1.0 


1600 


51 > 7 


HD 152248 


07Ib 


(48) 


RCW 106 


4681 -2375/+4871 


45.7 ± 10.0 


34.1 -8.1/+10.9 


2.5 


1100 


41 > 8 




05.5-6..5V 


(49) 


NGC 6823 


4983 -2584/+5685 


37.7 ± 4.0 


27.2 -5.2/+8.8 


2-4 


1900 


42 > 8 


BD 22° 3782 


07V(f) 


(50) 


RCW 121 


5323 -2671/+5390 


42.0 ± 4.0 


38.3 -8.3/+11.7 


4.2 


1600 


96 > 5 




06V-05V 


(51) 


NGC 2244 


5946 -3029/+6102 


68.9 ± 14.0 


54.1 -9.1/+10.9 


1.9 


1500 


54 (6 - 12) 


HD 46223 


04V(f) 


(52) 


NGC 2122 


6764 -3416/+6960 


74.7 ± 10.0 


58.8 -13.8/+11.2 


3.0 


48500 


52 > 9 


HD 270145 


061(f) 


(53) 


[OBS 2003] 179 


7000 ± 3500 


50.0 ± 20.0 


70.0 ± 10.0 


2-5 


7900 


10 > 16 


Obj 4 


Ofpe/WN9 


(54) 


Westerlund 2 


8845 -4456/+9009 


82.7" ± 5.5 


121.0 -43.8/+29.0 


1-2 


8000 


29 > 16.5 


WR20a A 


WN6ha 


(55) 



" Dynamical mass estimates from binary orbits exist for these stars. 



Designation -^ccl ''^max obs old "^max obs new D of StarS Id Wmcix Sp Type Rcf. 



[Mq] [Mq] [Mq] [Myr] [pc] 



RCW 95 


9670 


-4840/+9720 


87.6 ± 14.0 


64.6 


-9.6/+15.4 


1.5 


2400 


136 > 6 


- 


03V 


(56) 


IC 1805 


10885 


-5528/+11137 


75.6 ± 10.0 


56.5 


-6.5/+8.5 


2.0 


2350 


99 (6 - 12) 


HD 15558 


04-5in(f) 


(47) 


NGC 6357 


11978 


-6430/+11979 


115.9 ± 30.0 


79.8 - 


■14.8/+15.2 


1.0 


2560 


38 > 16.5 


HDE 319718A 


03If 


(57) 


NGC 3603 


1.3 


■ 10-* ± 3000 


150.0 ± 50.0 


121.0 


-41.0/+29.0 


0.7 


6000 




NGC 3603-B 


WN6ha 


(58) 


Trumpler 14/16 


17890 


-8945/+18676 


150.0 ± 50.0 


99.8 - 


■39.8/+50.2 


1.7 


2500 


64 > 16 


rj Carina 


LBV 


(59) 


NGC 6611 


25310 


-12659/+25503 


85.0 ± 15.0 


60.4 


-5.6/+9.6 


1.3 


1800 


460 > 5 


HD 168076 


04III 


(60) 


Cyg OB2 


75890 


-38453/+78716 


92.0 -25.0/+58.0 






2.0 


1700 


8600 > 1.3 


Cyg OB2-12 


B8Ia 


(61) 


Arches 


77225 


-39250/+77225 


120.0 -0.0/+30.0 


111.0 


-41.0/+39.0 


2.5 


7620 


196 > 20.1 


N4 


WN7-8h 


(62) 


R 136 


222912 - 


112104/+224426 


145.0 ± 15.0 


125.4 


-45.4/+24.6 


1-2 


48500 


8000 > 3 


R136al 


02If*/WN4.5 


(63) 



" Dynamical mass estimates from binary orbits exist for these stars. 
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APPENDIX C: NOTES ON INDIVIDUAL STAR CLUSTERS 
CI MWC 297 

Depending on t he uncertain distance the spectral type of the brightest star is either 09e (450 pc, Hillenbrand et al. 19921 ) or 



B1.5V (250 pc. lDrew et al.lll997l ). IWang fc LoonevI (|2007f ) assu me the 09e spectral type an d assign the star a mass of 26.5 



This mass e stimate is rather odd as older literature gives 22 ( Hanson et al.lll997l ) or 25.4 (jVacca et al.lll996l ) Mq for a 09 



star and the new iMartins et all |2005h calibration yields 18 Mq. Here we assume the newer distance of 250 pc and therefore 




Figure Bl. The yalues for the O star masses from the old calibration on the abscissae ys the newly calibrated ones on the ordinate. 
The thick solid line indicates were both yalues would be the same. 
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the spectral type B1.5V but include the mass of an 09V star as an upper limit (22 Mq for the old calibration and 18 Mq for 
the new one). 



C2 Pismis 11 

It is not certain whether HD 80077 is part of this star cluster or not ( Marco fc Neeueruela 20091 ) . If it is it would be one of 
the brightest stars in the MW. As it is an evolved star (B2Ia) the precise initial mass in not known, only that it should be 
above 40 Mq . 



C3 NGC 2244 

Bonatto &l Bica (|2009l ) give a mass of 625 Mq for this cluster from fitting a King pr ofile to the MS and PMS stars they found. 
But when using the number of BO to B3 stars (54 between 6 and 12 Mq Mq) from I Wolff et al.l (|2007l ). an IMF extrapolation 
yields a cluster mass of about 6000 Mq. A reason for this discrepancy is currently not known. 



C4 7 Velorum Cluster 

Jeffries et all (|2009l ) studied a group of several hundred PMS stars around the massive binary 7^ Vel (WC8 + 07.5V) in the 
Vela 0B2 association. They concluded that it is a young (~ 7 Myr) low-mass (250 - 360 Mq) cluster with a very massive star 
(initial mass for the WC8 star > 40 Mq). This object has not been included because of its age and the possible mass loss 
from the cluster due to gas expulsion. 



C5 Trumpler 14/16 

Trumpler 14 and T rumpler 16 are believed to be one single cluster optically divided by a large patch of gas and dust 
l Oev fc Claris I2OO5I ). Some studies do not consider rj Carina as part of this cluster. If 77 Car is not part of the cluster then 
the OSIf* star HD 93129AB (67^ Mq) would be the most-massive star in the cluster. 



C6 Cyg OB2 

The star Cyg OB2-12 is an evolved star (B8Ia) and is considered to be one of or maybe even the brightest star in the MW. 
Neeueruela et al. ( 20081 ) estimate its initial mass with 92 Mq. The region co ntains 8600 stars earlier than F3, some 2600 OB 
stars and around 120 O stars and has been described bv iKnodlsede ir i|2000l ) as a possible young globular cluster within the 
MW disk. 



C7 R 136 

There is some debate in the literature whether or not extremely ma ssive stars (> 8 Mq) might or might not look like 
Wolf-Rayet stars even while they still burn hydrogen in their cores (|Crowtherl 120071 ) . Therefore, the spectral type of the 
most-massive star (R136al) might be either 02If* or WN4.5. But this has little consequences for the mass estimate for this 
star as it is based on a detailed spectral analysis ( Massev fc Hunter 19981 ). The mass estimate for the cluster, 2.2 • 10^ Mq, 
is for R 136 only. The cluster contains about 8000 stars b etween 3 and 120 Mq and 39 03ir stars. The whole region of 30 
Doradus is believed to hold a total of 4.5 • 10^ Af© in stars ( Bosch et al. 20091 ). 



C8 CI 1806-20 

The cluster and its most massive star, LBV 1806-20, one of the brightest stars in the MW, have not been included because the 
cluster is known to host a pulsar (SCR 1806-20) and is therefore most probably too old for the current study. Additionally, 
very little is known about its stellar population besides the LBV, the pulsar, an OB sup ergiant and three Wolf-Rayet stars 
as this cluster is on the far side of the Calaxy (D 15000 pc) and heavily dust obscured (|Figer et al.ll2005l ). 



C9 Quintuplet 

This massive cluster (~ 1.2-10'' Mr^. lPiger et al.ll 19991 : [Liermann et al.ll2009l ) hosts the Pistol star which is one of the brightest 
and most massive stars in the Galaxy (150^ Mq?) but is believed to be 4-6 Myr old and therefore too old for this study. 



22 C. Weidner, P. Kroupa and I. A. D. Bonnell 



CIO Westerlund 1 



Another massive cluster (> 10* Mq, Brandner 20081 ') with about 150 O stars. It houses an X-ray pulsar I Muno et al. 20061 ) 



and is therefore not included into this study as it is too old (4-5 Myr) 
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